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Abstract. Deep 21-cm Hi line observations of the Virgo cluster spiral galaxy NGC 4569 have been obtained with 
the VLA in its D configuration and with the Efi'elsberg 100-m telescope. A low surface density arm was discovered 
in the west of the galaxy, whose velocity field is distinct from that of the overall disk rotation. The observed 
gas distribution, velocity field, and velocity dispersion are compared to snapshots of dynamical simulations that 
include the effects of ram pressure. Two different scenarios were explored: (i) ongoing stripping and (ii) a major 
stripping event that took place about 300 Myr ago. It is concluded that only the post -stripping scenario can 
reproduce the main observed characteristics of NGC 4569. It is not possible to determine if the gas disk of 
NGC 4569 had already been truncated before it underwent the ram pressure event that lead to its observed Hi 
deficiency. 
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1. Introduction 

The spiral galaxy NGC 4569 is a very peculiar member of 
the Virgo cluster. Its very large size {D25—9'.5, or 47 kpc 
at the assumed Virgo cluster distance of 17 Mpc) com- 
pared to other cluster members, and negative radial ve- 
locity {vi- = —235 kms^^ with respect to the cluster mean 
velocity of ^ 1100 kms^^) have lead to some doubt about 
its cluster membership, despite its close projected distance 
to the cluster center {d = 1.7°). Van den Bergh (1976) 
classified this very red galaxy as anemic due to its low 
arm~interarm contrast on optical images. H i line synthesis 
observations (Warmels 1988; Cayatte et al. 1990) showed 
that it has lost more than 90% of its H i gas and that its 
H I distribution is heavily truncated. NGC 4569 is thus an 
exceptional galaxy, it is one of the largest and most H i de- 
ficient galaxies in the Virgo cluster (Solanes et al. 2001). 
Since field galaxies are generally not Hi deficient, this is 
the proof that NGC 4569 is a Virgo cluster member whose 
atomic gas has been most probably stripped due to its 
rapid motion through the hot intracluster medium (ICM) . 
Nevertheless, it is not clear when and where NGC 4569 
lost its gas. 
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A possible hint on its stripping history may be derived 
from X-ray observations. The analysis of ROSAT obser- 
vations (Tschoke et al. 2001) showed that there is no soft 
emission (0.1-0.4 keV) from the northern half of the disk, 
whereas it is pronounced in the southern disk. This is con- 
sistent with a scenario where NGC 4569 is moving to the 
north-east, which causes ram pressure to push the hot ISM 
in the galaxy towards the south-west. 

In addition, NGC 4569 harbors a bright, compact op- 
tical nucleus that , according to Keel (1996), is due to a 
starburst dominated by supergiant stars. Tschoke et al. 
(2001) suggest that a part of the observed extended X-ray 
emission is due to an outflow of hot ISM from the nu- 
cleus. The presence of such a feature in the west and its 
absence in the east is also consistent with a galaxy mo- 
tion towards the north-east. In the following we will not 
discuss the possiblity of nuclear outflow, since it expands 
vertically to the disk and does not affect the dynamics of 
the atomic gas, which is mainly located in the disk. 

In the present paper we report on new deep VLA 21 cm 
Hi observations of NGC 4569, which reveal a new and 
unexpected feature - a kinetically perturbed, low surface 
brightness arm - and we present models of the gas distri- 
bution and velocity field based on numerical simulations 
that include the effect of ram pressure stripping. This en- 
ables us to date the latest ISM-ICM interaction that has 
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Table 1. Physical parameters of NGC 4569 



Other names M 90 

VCC 1690 
UGC 7786 

a (2000)'' 12''36'"49.8' 

& (2000)^ +13°09'46" 

Morphological type^ Sab 

Distance to the cluster center (°) 1.7 

Optical diameter D25'' (') 9.5 

B?,'' 9.79 

Systemic heliocentric velocity'' (kms~^) — 235±3 

Distance D (Mpc) 17 

Vrotmax (kms-i) 250*^ 

PA 24° ^ 

Inclination angle 56° '',62° 

HI deficiency^ 1.3±0.2 



RC3, de Vaucouleurs et al. (1991) 
^ based on Hi inematics, Guharthakurta et al. (1988) 
° Cayatte et al. (1994) 

based on optical photometry, Bottinelli et al. (1983) 



led to the large Hi deficiency of NGC 4569. The present 
work is a part of a series of articles where we aim to un- 
derstand the physics of ICM-ISM interactions (VoUmer 
et al. 1999, VoUmer et al. 2001, VoUmer 2003, VoUmer & 
Huchtmeier 2003). 

This article is organized in the following way: the ob- 
servations and their results are described in Sect. |2] and 
Sect. 13 respectively, the dynamical model is presented in 
Sect. 01 our model snapshots are compared with the obser- 
vations in Sect. Inland discussed in Sect.|^ The conclusions 
are given in Sect.[7| 

2. Observations 

2.1. VLA 

The 21 cm line observations were made with the Very 
Large Array (VLA), for description see Napier et al. 
(1983). The field was centered on NGC 4569 (see Tab.P). 
We observed on November 23 and 29 2001 for a total of 
290 minutes with the D configuration. The observed band- 
width of 3 MHz was devided into 63 velocity channels of 
48.8 kHz (~ 10 kms~^) each. The velocity channels are 
centered on Whoi ~ —235 kms^^. The data were calibrated 
using the standard VLA data reduction program package 
(AIPS). All strong continuum sources were subtracted di- 
rectly in the UV plane using a linear interpolation of the 
UV data points of the first and last 10 channels. The fi- 
nal rms noise per 10 kms~^ wide channel at a smoothed 
spatial resolution of 52" x 40" is 0.6 mJy/beam, which is 
close to the expected value. We combined the UV data 
directly with the C array data of Cayatte et al. (1990), 
which were processed in the same way. The same weight 
was attributed to both data sets. Since the C array data 



Table 2. Integration times and rms. 



position 


C 


NW 


W 


SW 


SE 


NE 


At (min) 


120 


120 


120 


120 


120 


120 


rms (mjy) 


1.3 


1.1 


2.0 


1.6 


1.6 


1.5 



have a channel width of 20 kms""', we binned every 2 
channels of our D array data before joining it with the C 
array data. The velocity offset between the two data sets is 
Av — 4 km s^^ . The resulting image of the joined UV data 
was CLEANed with a 39" x 31" FWHM beam. We ended 
up with a r.m.s. noise of 0.5 mJy/beam in one 20 kms^^ 
wide channel. A flux of 1 mJykms^^/beam corresponds 
to an Hi column density of 10^* cm^'^. 

2.2. Effelsberg 

On March 14-21, 2002, we performed 21-cm line obser- 
vations at 6 different positions centered on the systemic 
velocity of NGC 4569 with a bandwidth of 12.5 MHz. 
The two-channel receiver had a system noise of ~30 K. 
The 1024 channel autocorrelator was split into four banks 
with 256 channels each, yielding a channel separation of 
~10 kms^^. We further binned the channels to obtain 
a final channel separation of ~20 kms~^ like the VLA 
data. The galaxy's central position and four positions at 
a distance of one beam width (9.3') to the NW, SW, SE, 
and NE from the galaxy center were observed in on-off 
mode (5 min on source, 5 min off source). In addition, we 
observed a sixth position 6^5 west of the galaxy center. 
Care was taken to avoid other Virgo cluster galaxies with 
velocities within our bandwidth in all observations. We 
used 3C286 for pointing and flux calibration. The obser- 
vation time was 120 min per position. The resulting noise 
(Table is partly determined by small amplitude inter- 
ferences, but it is close to the theoretical noise of 2 mJy 
per hour of integration: on average 1.5 mJy (varying from 
1.1 to 2.0 mJy). 

3. Results 

3.1. The VLA data 

The channel maps are presented in Fig. ^ Although two 
channels {v — —30 kms^^ and —10 kms^^) are contam- 
inated by Galactic Hi, the Hi of NGC 4569 at these ve- 
locities can still be recognized and isolated from it. The 
galactic rotation of NGC 4569 is clearly visible. The most 
interesting feature appears in the velocity range between 
—400 kms""' and —276 kms^^, where a detached emission 
region can be seen in the south-west part of NGC 4569, 
which shows rotation. 

We have produced an integrated H i spectrum of 
NGC 4569 by adding all signals exceeding the 3(t level 
in all channels that contain line emission, inside an area 
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Fig. 1. VLA Hi channel maps of NGC 4569. The central velocity of each channel is indicated in the upper part of 
each panel. The HPBW (31" x 39") is plotted in the lower left corner of the first channel. The contours are in steps 
of 3(7=1.5 mjy/beam beginning at 1.5 mJy/beam. 



around NGC 4569 small enough not to pick up Galactic 
H I. Fig. 13 shows the integrated VLA spectrum together 
with the Effelsberg 100-m data of the central pointing. The 
Effelsberg beam of 9.3' includes the the whole disk Hi 
emission. The Galactic foreground is responsible for the 
huge peak in the Effelsberg spectrum at —220 kms^^ 



relative to the systemic velocity, whose width (FWHM) 
is Aw ^ 50 kms~^. The total flux and peak flux density 
of Cayatte et al. (1990) are --20% higher than our VLA 
(C and D array) values, which is within the calibration 
accuracy. Both of our spectra agree reasonably well, and 
both peaks of the double horned shaped spectrum have 



4 



VoUmer et al.: NGC 4569 



NGC 4569 
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Fig. 2. Solid line: EfFelsberg spectrum of the central po- 
sition. Dashed line: integrated VLA spectrum. Dotted 
line: 3a noise level of the Effelsberg spectrum. Plotted 
are radial velocities relative to the systemic velocity of 
NGC 4569. 



about the same peak flux density. We derive a linewidth 
at 20%/50% of the peak flux of W2Q = 390 kms-i/ 
W50 = 350 kms~^ and a total flux of 8.85 Jykms"^, 
which corresponds to a total H i mass of 6 10^ Mq assum- 
ing a distance of 17 Mpc. 

The first three moment maps (column density, velocity 
field, and dispersion velocity) can be seen in Fig. 13 For the 
determination of all moment maps the data were clipped 
at a level of 3a = 1.5 mJy/beam in each channel con- 
taining Hi emission. The Hi distribution consists mainly 
of three maxima: at the center, to the north and to the 
south. The southern maximum lies west from the major 
axis. The H i is thus not uniformly distributed within the 
disk, but might be distributed along spiral arms as it is 
the case for NGC 4579 (Cayatte et al. 1990) or NGC 4548 
(VoUmer et al. 1999). We observe a low surface density 
arm to the west, which was already visible on the channel 
maps. This arm can also be seen on the Ha image pre- 
sented in Tschoke et al. (2001) and the blue optical image 
shown in Fig. 0] 

The velocity field of NGC 4569 is quite symmetric 
within the disk. The roation curve is rising on both sides 
with a maximum of 250 kms~^ assuming the kinematic 
inclination angle oii = 56° (Guharthakurta et al. 1988) or 
230 kms^^ assuming the inclination based on the optical 
morphology i = 62° (Bottinelli et al. 1983). The western, 
low surface density arm is kinematically distinct from the 
disk gas. This arm shows rotation, but has systematically 
higher velocities than the disk gas. 

As expected, the velocity dispersion shows a maximum 
in the galaxy center, due to rotation and beam smearing. 
The absolute maximum is located ~ 10" to the north-east 
of the nucleus and has a velocity dispersion greater than 
80 kms""'^, a second more extended maximum is observed 
at the inner edge of the western, low surface density arm 
south of the nucleus. 




RIGHT ASCENSION (B195D) 




Fig. 3. Hi moment maps of NGC 4569. Upper panel: 
column density. The contour steps are from 1 to 9 x 
1.3 10^° cm~^. Middle panel: velocity field. Contour levels 
are from -440 to -40 kms^^ in steps of 20 kms~^. Lower 
panel: dispersion velocity. Contour levels are from 1 to 
9 X 8.3 kms-^ 



3.2. The EfFelsberg data 

For the 4 off-center Effelsberg positions we have synthe- 
sized spectra from the VLA data using the area defined 



VoUmer et al.: NGC 4569 



5 



13 30 
29 
28 

o 
m 

S 27 

z 
O 

i 26 

O 
LU 
Q 

25 




_l I \ \ L 

25 20 15 10 05 

RIGHT ASCENSION (B1950) 



Fig. 4. Contours: Hi column density. The contour levels 
are the same as in Fig. 01 Greyscale: B band DSS image. 
The contrast is chosen such that the whole disk extent 
becomes visible. 

by the pointing positions and the (9f3) EfFelsberg HPBW. 
In this way we take advantage of the spatial H i distribu- 
tion of the VLA data (see Vollmer & Huchtmeier 2003). 
Fig.Elshows the Effelsberg spectra of these positions (solid 
lines) together with the synthesized VLA spectra. Since 
the VLA data cube shows only H i gas located in the disk 
(see Cayatte et al. 1990), any emission excess in the off- 
center Effelsberg spectra must be due to gas located out- 
side the disk. No such excess is seen and an upper mass 
limit ~ 10^ M0 can be set for the Hi outside the disk, 
corresponding to a signal of 5a in a 20 kms~^ channel 
with the an rms of a ~ 1.6 niJy. The absence of atomic 
gas far away from the galaxy (> 20 kpc), as was found for 
NGC 4388 (Vollmer & Huchtmeier 2003), is an important 
constraint on dynamic stripping models. 

4. Simulations 

We have adopted a model where the ISM is simulated 
as a collisional component, i.e. as discrete particles which 
possess a mass and a radius and which can have inelas- 
tic collisions (sticky paticles). Since the ISM is a turbu- 
lent and fractal medium (see e.g. Elmegreen & Falgarone 
1996), it is neither continuous nor discrete. The volume 
filling factor of the different phases is smaller than one. 
The warm neutral and ionized gas fill about 50% of the 
volume, whereas cold neutral gas has a volume filling fac- 
tor of only 2% (Boulares & Cox 1990). It is not clear how 
this fraction changes, when an external pressure is ap- 
plied. In contrast to smoothed particles hydrodynamics 
(SPH), which is a quasi continuous approach and where 
the particles cannot penetrate each other, our approach 





Fig. 5. Solid lines: Effelsberg 100-m spectra of the five 
off-center positions. Their locations with respect to the 
galaxy center are marked on top of each panel. Dashed 
line: synthesized VLA spectra, which only show Hi disk 
emission. Dotted line: 3a noise levels of the Effelsberg 
spectra. Radial velocities are given relative to the systemic 
velocity of NGC 4569. 



allows a finite penetration length, which is given by the 
mass-radius relation of the particles. Both methods have 
their advantages and their limits. The advantage of our 
approach is that ram pressure can be included easily as 
an additional acceleration on particles that are not pro- 
tected by other particles (see Vollmer et al. 2001). In this 
way we avoid the problem of treating the huge density 
contrast between the ICM (n ^ lO"^ cm"^) and the ISM 
(n > 1 cm"'^) of the galaxy. 

Since the model is described in detail in Vollmer et al. 
(2001), we will summarize only its main features. The N- 
body code consists of two components: a non-coUisional 
component that simulates the stellar bulge/disk and the 
dark halo, and a collisional component that simulates the 
ISM. The 20 000 particles of the collisional component rep- 
resent gas cloud complexes which are evolving in the grav- 
itational potential of the galaxy. 

The total assumed gas mass is Mg°* = 6.210^ Mq, 
which corresponds to the total neutral gas mass before 
stripping, i.e. to an Hi deficiency of 1.0, which is defined 
as the logarithm of the ratio between the Hi content of 
a field galaxy of same morphological type and diameter 
and the observed Hi mass. To each particle a radius is 
attributed depending on its mass. During the disk evolu- 
tion the particles can have inelastic collisions, the outcome 
of which (coalescence, mass exchange, or fragmentation) 
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Table 3. Total mass, number of particles N, particle mass 
M, and smoothing length I for the different galactic com- 
ponents. 



component 


Mtot (Mq) 


N 


M (Mo) 


I (pc) 


halo 


2.4 10^' 


16384 


1.410' 


1200 


bulge 


8.210^ 


16384 


5.0 lO'' 


180 


disk 


4.1 10^" 


16384 


2.510*^ 


240 



is simplified following Wiegel (1994). This results in an 
effective gas viscosity in the disk. 

As the galaxy moves through the ICM, its clouds are 
accelerated by ram pressure. Within the galaxy's inertial 
system its clouds are exposed to a wind coming from a 
direction opposite to that of the galaxy's motion through 
the ICM. The temporal ram pressure profile has the form 
of a Lorentzian, which is realistic for galaxies on highly 
eccentric orbits within the Virgo cluster (VoUmer et al. 
2001). The effect of ram pressure on the clouds is sim- 
ulated by an additional force on the clouds in the wind 
direction. Only clouds which are not protected by other 
clouds against the wind are affected. 

The non-collisional component consists of 49 125 par- 
ticles, which simulate the galactic halo, bulge, and disk. 
The characteristics of the different galactic components 
are shown in Table |2I The resulting rotation velocity is 
~220 kms~^ and the rotation curve is flat. 

The particle trajectories are integrated using an adap- 
tive timestep for each particle. This method is described 
in Springel et al. (2001). The following criterion for an 
individual timestep is applied: 



20 kms~^ 



(1) 



where Oi is the acceleration of the particle i. The mini- 
mum value of ti defines the global timestep used for the 
Burlisch-Stoer integrator that integrates the coUisional 
component. 

The galaxy is on an eccentric orbit within the cluster. 
The temporal ram pressure profile can be described by: 



Pn 



f2 -I- ^2 



(2) 



where inw is the width of the profile (VoUmer et al. 2001). 
We set pmax=4000 cm-3(km/s)2 and tew=50 Myr. The 
efficiency of ram pressure also depends on the inclination 
angle i between the galactic disk and the orbital plane 
(VoUmer et al. 2001). We set i=35°. The model parameters 
are listed in Tabled 

The large scale evolution of NGC 4569 can be seen in 
Fig. El As expected, the stellar disk does not change, be- 
cause ram pressure selectively affects only the gas. During 
the passage of NGC 4569 through the cluster core {t = 
Myr corresponds to the minimum distance between the 
galaxy and M87) its gas is pushed strongly towards the 



Table 4. Model parameters for time dependent ram pres- 
sure stripping 



maximum ram pressure (cm~^(kms~^)^) 4000 
inclination angle between orbital and disk plane 35° 
final HI deficiency 1.0 



face— on 



PA=23' 



65" 




-20 -10 10 20 




-20 -10 10 20 
RA (kpc) 



Fig. 7. Comparison snapshot of the gas distribution of 
NGC 4569 for the case of ongoing ram pressure stripping. 
Left panel: the disk seen face-on. Right panel: view of the 
disk using the observed position angle and inclination of 
NGC 4569. The arrow indicates the direction of the wind, 
i.e. opposite to the galaxy's motion. 



south-west and finally stripped at i > 50 Myr. The gas 
that was pushed away, but not accelerated to the escape 
velocity, falls back onto the galaxy at i > 150 Myr. This 
re-accretion of the gas has to be regarded with caution, 
since its occurence depends on the gas physics (cloud ex- 
pansion, evaporation, ISM-ICM phase mixing) used in the 
model. We assume that the column density of the clouds 
does not change during the interaction - however, if it 
were to decrease due to, e.g., evaporation or ICM-ISM 
mixing, then re-accretion will be suppressed. 

During the whole galaxy evolution within the cluster 
when ram pressure is important only two scenarios can 
give rise to an asymmetric gas arm: (i) during the reset- 
tlement of the gas after a strong stripping event and (ii) 
during the phase of ongoing stripping. Therefore, we com- 
pare two snapshots with our Hi observations: (i) the last 
snapshot of Fig. El (i = 300 Myr), i.e. a past stripping 
event and (ii) the case of ongoing stripping - another sim- 
ulation made using a truncated gas disk at i? = 10 kpc. 
withp,nax=2000cm-3(km/s)2, i ^ 35°, iHW==80Myr. The 
comparison snapshot at i = Myr is shown in Fig. [71 For 
a given position angle and inclination, the third angle, 
which is the azimuthal angle in the plane of the galaxy, 
has to be fixed. A variation of the azimuthal angle leads 
to a variation of the 3D direction of the galaxy's motion 
with respect to the observer. We chose the value that re- 
produces best the observed, western, low surface density 
arm. In both cases the galaxy moves to the north-east. 
We made sure that the galaxy's radial velocity did not 
faU below -1500kms-i. 
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Fig. 6. Evolution of the model stellar (1st and 3rd column) and gas disk (2nd and 4th column). The major axis 
position angle and inclination of NGC 4569 are PA = 23° and i = 62°, respectively. The arrow indicates the direction 
of ram pressure, i.e. it is opposite to the galaxy's velocity vector, and its size is proportional to pw|a,i- The galaxy 
passes the cluster core at Myr. The time to the core passage is marked at the top of each panel. 



We constructed model data cubes of the two snapshots 
and made moment maps with AIPS. Only gas particles 
with a local gas density greater than 100 cm~'^ where put 
into the model cube. 

5. Comparison between simulations and 
observations 

In this Section we compare model spectra and the three 
model moment maps (gas distribution, velocity field, and 
dispersion velocity) to those of our H i observations (Fig.|21 



13 andlSj. The initial velocity field of the gas disk is shown 
in Fig. 121 To facilitate the comparison with the final snap- 
shots of our simulation, we have truncated the disk at a 
radius of 5 kpc. Our model velocity field has a steeper ro- 
tation curve and reaches the plateau of constant rotation 
velocity at a smaller radius than the observed one (Fig.lSJ. 
However, the CO position- velocity diagram along the ma- 
jor axis of NGC 4569 (Fig. 23 of Kenney & Young 1988) 
shows a steeper rotation curve, which is close to that of 
our model. 
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Fig. 8. Initial gas velocity field. The disk is truncated at 
a radius of 5 kpc. 
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Fig. 9. Model H i spectrum for ongoing ram pressure strip- 
ping. 

5.1. Ongoing stripping 

The model spectra corresponding to the observed cen- 
ter spectrum and those of the four offset H i observations 
(Fig.|5l) for ongoing stripping are shown in Fig. l^andlTUl 
The two peaks of the central double-horn spectrum have 
about the same velocity width and the edges of the profile 
are sharp. The off-center spectra in the NE, SE, and NW 
show maximum flux densities smaller than the detection 
limit of our Effelsberg data (~ 5 mJy). The signature of 
the western arm is clearly visible at negative velocities in 
the W and SW spectra. 

Fig. llll shows the moment maps for ongoing ram pres- 
sure stripping. As already mentioned, the projection pa- 
rameters are chosen such that (i) the radial velocity of 
NGC 4569 and (ii) the observed, western, low surface den- 
sity arm are reproduced. The central hole in the gas dis- 
tribution is artificial, in order to avoid too small timesteps 
for the particle integration. The north-eastern maximum 
represents compressed gas that has been pushed by ram 
pressure to smaller galactic radii. The velocity field in the 
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Fig. 10. Model spectra corresponding to our offset Hi ob- 
servations (Fig. ISJ for ongoing stripping. 

north is that of a rising rotation curve, because the ISM 
is pushed inwards by the ICM. This effect is also observed 
for NGC 4654 (VoUmer 2003). The southern part of the 
disk velocity field shows a plateau of constant rotation ve- 
locity. The velocity contours are continuous from the disk 
to the western, low surface density arm, because the gas 
there is accelerated and pushed out of the galaxy in the 
sense of rotation. As expected, the maximum of the ve- 
locity dispersion is located in the galaxy center. A second 
maximum can be found in the north-west. It corresponds 
to gas that has been pushed inwards in the north-east by 
the ICM and has rotated (the galaxy is rotating clock- 
wise) towards the west. 

5.2. Past stripping event 

The model spectra corresponding to the observed cen- 
ter spectrum and those of the four offset H i observations 
(Fig. O for a past stripping event are shown in Fig. E| 
and^l The low velocity peak of the central double- horn 
spectrum has a larger velocity width than the high velocity 
peak. Moreover, the low-velocity edge of the profile shows 
a pronounced wing. All off-center spectra except the one 
in the W show maximum flux densities smaller than the 
detection limit of our Effelsberg data (~ 5 mJy). 

Fig. 1141 shows the moment maps for a past ram pres- 
sure stripping event. Again, the projection parameters are 
chosen such that they reproduce (i) the radial velocity 
of NGC 4569 and (ii) the observed, western, low surface 
density arm. Again, the central hole in the gas distribu- 
tion is artificial, in order to avoid too small timesteps for 
the particle integration. The spatial extent of the model 
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Fig. 12. Model H i spectrum for a past ram pressure strip- 
ping event. 
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Fig. 11. Model moment maps for ongoing ram pressure 
stripping. Upper panel: gas distribution in arbitrary units. 
The relative contours are the same as in Fig. |3| Middle 
panel: velocity field. The contours are the same as in 
Fig. 121 Lower panel: velocity dispersion. The contours are 
the same as in Fig. El 




Fig. 13. Model spectra corresponding to our offset Hi ob- 
servations (Fig. 01 for a past ram pressure stripping event. 

gas distribution is smaller than that of the observed one, 
because we have decided to reproduce the Hi deficiency 
rather than the spatial extent. For simplicity we assume 
that all gas in the model is atomic. Consequently, our ini- 
tal central gas surface density is higher than the observed 
one (our model rather represents an Sc/Sb galaxy, whereas 
NGC 4569 is classified as Sab). Thus, for a given Hi de- 
ficiency, our model gas disk is smaller than the observed 
disk. The disk velocity field shows a plateau on both sides, 
which is reached at smaller galactic radii in the south than 
in the north. The velocity field of the southern part of the 
western, low surface density arm is not a continuation of 
that of the disk. The absolute velocities of the arm are 
higher than those of the neighbouring disk. The maximum 
of the velocity dispersion is located in the center. A sec- 
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ond, point-like maximum can be found to the south-east 
of the nucleus where the radial velocity has its minimum 
value. This feature is due to a gas streamer that falls onto 
the galaxy from behind. In addition, we observe an ex- 
tended maximum at the inner, edge of the western, low 
surface density arm south of the nucleus. 

6. Discussion 

Since both models reproduce the observed, western, low 
surface density arm as well as the radial velocity of 
NGC 4569, it became clear that it is impossible to chose 
between the two scenarios only on the basis of the ob- 
served gas distribution. Only the velocity information in 
addition to the gas distribution allows us to discriminate. 

The extended low velocity edge of the observed cen- 
tral H I spectrum is only reproduced by the past stripping 
simulation. In general, the offset spectra of both simula- 
tions are consistent with our observations, i.e. only in the 
west (W) atomic gas is detected (flux density greater than 
~ 5 mjy in one 20 kms~^ channel). The western arm of 
the ongoing stripping model is somewhat too prominent 
in the model SW spectrum, but this is due to the detailed 
choice of the stripping parameters. It is important to note 
that the tail of stripped and back-falling gas of the post- 
stripping model (Fig.ini) could not have been detected with 
the Effelberg telescope (Fig. neither with the VLA). 

Ongoing stripping steepens the velocity field in the di- 
rection of the galaxy's motion and elongates the velocity 
contours on the opposite side (Fig. ^2 see also VoUmer 
2003), because ram pressure pushes the gas to smaller 
galactic radii in the direction of the galaxy's motion and 
accelerates it in the sense of rotation on the other side. 
On the other hand, a past stripping event causes the 
steepening and asymmetry of the rotation curve (Fig. 1141) . 
This asymmetry is also seen in the observed velocity field 
(Fig. EJ where the rotation curve of the souther half is 
steeper than that in the north. Regarding the general 
steepening of the rotation curve: since we do not know 
the initial velocity field of NGC 4569, we do not think 
that the difference between the simulated (Fig. I14|) and 
observed (Fig. velocity fields poses a problem for the 
validity of the simulations. 

Whereas the velocity field in the western arm region 
is continuous in the case of ongoing stripping, it shows 
a discontinuity in the post-stripping scenario as is ob- 
served. Also, the difference in absolute velocity between 
the disk and the western arm of the post-stripping sce- 
nario fits nicely our H i observations. As for the distri- 
bution of dispersion velocities, the observed, single, local, 
extended maximum at the inner edge of the western am 
is only reproduced by the post-stripping model. It is lo- 
cated at the interface between the western, low surface 
density arm and the gas rotating within the galactic disk. 
We failed to find any suspect features in the 3D gas dis- 
tribution at that location. The general dynamics are very 
complicated, which makes it hard to disentangle different 
physical effects due re-accretion, gas expansion, and rota- 



tion. The lines in this region are broadened to both larger 
and smaller radial velocities compared to a region south 
of it. We have no clear explanation why this region shows 
an enhanced dispersion velocity. 

There is also a point-like maximum visible, like the 
one observed, but it is located to the south. The enhanced 
linewidth arises from a back-falling gas streamer. Since the 
exact dynamics of such a streamer depends on gas physics 
(e.g. volume filling factor, shock heating, evaporation) and 
the detailed stripping parameters, we do not expect to 
fit its exact location. We can only speculate that such a 
region of unusual large velocity dispersion might be a sign 
of re-accretion. 

We thus conclude that the ongoing stripping scenario 
is ruled out based on the comparison of the observed and 
simulated velocity fields and velocity dispersion distribu- 
tions. On the other hand, the post-stripping scenario re- 
produces the major characteristics of the observations: the 
truncated gas disk, the western, low surface density arm, 
and the discontinuity between the velocity field in the 
galactic disc and the western arm. In order to investigate 
if the gas disk of NGC 4569 was already truncated before 
it underwent the ram pressure stripping event that lead 
to the observed H i deficiency, we made an additional sim- 
ulation starting with a truncated gas disk at i? = 10 kpc. 
The resulting moment maps are not distinguishable from 
that of Fig. El Thus we conclude that it is not excluded 
that NGC 4569 has already experienced a ram pressure 
event more than ~ 1 Gyr ago that lead to a truncated 
H I disk. After such a first ram pressure stripping event 
the gas disk had enough time to evolve into an equilib- 
rium state without major perturbations. Possibly, the gas 
disk became more robust via the 'disk annealing process' 
described by Schulz & Struck (2001). This is consistent 
with our finding that the rotation curve steepens after the 
ram pressure event. During the galaxy's way through the 
outer parts of the cluster some of its ISM might have been 
replenished by the released gas by dying stars. All we can 
say is that NGC 4569 had a less truncated gas disk be- 
fore the last ram pressure stripping event than it has now. 
For each stripping event, the following parameters are im- 
portant for the efficiency of ram pressure stripping: the 
gas disk truncation radius, the galaxy orbit, and the in- 
clination angle between the galaxy's disk and the orbital 
plane during the cluster core passage. Different core pas- 
sages can have different maximum ram pressures (Vollmer 
et al. 2001), and the one with the highest maximum ram 
pressure is decisive for the galaxy. NGC 4569 experienced 
its most efficient ram pressure stripping event ~ 300 Myr 
ago. 

A final comparison of the data cubes can be directly 
made in 3 dimensional representations. We have developed 
a 3D visualisation of data cubes which enables us to eas- 
ily distinguish different kinematic features (see Vollmer et 
al. 1999). All velocity channels seen in Fig. ^ are piled up 
to give the cube, whose axes are the right ascension, the 
declination, and the radial velocity. All points in the cube 
having intensities exceeding a chosen level (2 mJy/beam) 
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become opaque, the rest remain transparent. The surface 
created in this way is iUuminated by hght which is com- 
ing out of the observer's direction, making brighter fea- 
tures look closer to the observer. If we had, for example, 
a sphere of increasing intensity to the centre, it would ap- 
pear as a ball with its radius increasing with decreasing 
chosen intensity level. This representation allows to anal- 
yse the whole cube from any possible point of view. In 
Fig. El we have rotated the cube in a way that the RA- 
Dec plane is perpendicular to the image plane, and the 
velocity axis is the vertical axis. The kinematics of rota- 
tion form a disk-like structure in the restricted phase space 
(RA, DEC, Wr). We have rotated the cube such that this 
disk structure is seen edge-on. The observed, western, low 
surface density arm can be seen in the upper right part of 
the image and is clearly distinct from the galaxy rotation 
fFig. 115b ). For comparison we show the unperturbed gas 
disk in Fig. 115b . The western arm in the ongoing strip- 
ping scenario (Fig. 1151 ) is spatially more extended and 
extends to smaller velocities. That in the post-stripping 
scenario (Fig. 1151 11 nicely fits the 3 dimensional extent of 
the observed western arm. However, due to the southern 
plateau of the model velocity field, the southern part of the 
disk-like structure in the restricted phase space is thicker 
than the observed one. We thus conclude that the post- 
stripping scenario is the most probable one for NGC 4569. 

This galaxy is not the first one, where an extended, 
low surface density arm has been discovered in deep H i 
imaging observations, as an extended, low surface density 
arm has also been found in three other Virgo spiral galax- 
ies: NGC 4321 (Knapen et al. 1993), NGC 4654 (Phookun 
& Mundy 1995) and NGC 4548 (VoUmer et al. 1999). An 
extreme case is NGC 4388 where VoUmer & Huchtmeier 
(2003) detected neutral hydrogen more than 20 kpc above 
the disk plane. These low surface density features stand 
in contrast to high surface brightness, extraplanar emis- 
sion regions, which are located close to the disk. These 
have been observed in NGC 4522 (Kenney et al. 2003) 
and NGC 4438 (Cayatte et al. 1990), two galaxies that 
are most probably actively stripped. 

The question arises if the high and low surface density 
features are caused by two different stages of an ICM- 
ISM interaction. We might speculate that the low surface 
brightness features whose velocity fields are not a continu- 
ation of that of the neighbouring disk are reminiscent of a 
past stripping event (NGC 4569 and NGC 4548), whereas 
those with a continuous velocity field are due to weak 
ongoing ram pressure stripping, perhaps in addition to a 
gravitational interaction (NGC 4654). On the other hand, 
high surface brightness features appear to be evidence of 
a major ongoing stripping process. 

If we assume a temporal ram pressure profile as pre- 
sented in VoUmer et al. (2001) the observational groups 
might be understood in the following way: 

(i) strong stripping leading to a large Hi defi- 
ciency (NGC 4388, NGC 4438, NGC 4522, NGC 4548, 
NGC 4569): during the maximum ram pressure phase 



high surface density gas is pushed out of the galactic disk 
(NGC 4438, NGC 4522). When the pushed-out gas moves 
to larger distances from the galaxy it expands and evap- 
orates, leading to low surface density Hi (NGC 4388). 
During active ram pressure stripping high surface density 
gas is accelerated and pushed to smaller galactic radii in 
the direction of the galaxy's motion and to larger galactic 
radii at the opposite direction. When the galaxy leaves 
the cluster core, the ram pressure decreases strongly and 
the gas is left with just the gravitational potential of the 
galaxy, to which it tries to adapt itself. The compressed 
region expands and the pushcd-away gas falls back to the 
galaxy. This creates expanding, asymmetric ring struc- 
tures whose surface density decreases during their expan- 
sion. These short-lived (< 10^ yr) features represent the 
kinetically perturbed low surface density arms observed in 
NGC4569 (and maybe NGC 4321 and NGC4548, whose 
cases are less clear). The number of observed low surface 
density arms in Virgo cluster spirals is small, but since 
these observations are time consuming, only few galaxies 
have been observed deep enough to allow their detection. 
However, all observations of sufficient depth have shown 
the presence of a low surface density arm. 

(ii) weak stripping together with a gravitational in- 
teraction leading to a small or even negligible Hi defi- 
ciency: ram pressure pushes away the low surface density 
gas which has been torn out of the galactic disk by the 
gravitational interaction (NGC 4654). 

This hypothesis needs to be verified with a larger sam- 
ple of cluster spirals. 

7. Conclusions 

We obtained deep H i line observations of the anemic Virgo 
spiral galaxy NGC 4569 with the VLA in its D configura- 
tion and with the Effelsberg 100-m telescope. Snapshots 
of a dynamical model including the effects of ram pressure 
were compared to the observations. Two different stages 
of a ram pressure stripping event were investigated: (i) on- 
going stripping and (ii) a past stripping event. The results 
of our investigations are: 

1. we detected a low surface density arm in the west of 
the galaxy; 

2. the velocity field of this arm is not a continuation of 
the velocity field of disk rotation; 

3. the absolute radial velocities of the arm are higher than 
those of the neighbouring disk; 

4. the velocity dispersion shows a local maximum at the 
inner edge of the arm; 

5. there is no gas mass greater than 10'' Mq detected far 
away (> 20 kpc) from the galaxy; 

6. the post-stripping scenario can reproduce the main 
characteristics of the gas distribution, velocity field, 
and velocity dispersion together with the radial veloc- 
ity of the galaxy; 

7. within this scenario the galaxy has passed through the 
cluster core about 300 Myr ago; 
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8. it is not excluded that NGC 4569 had an ah'eady trun- 
cated gas disk before this ICM-ISM interaction. 
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Fig. 14. Model moment maps for a past ram pressure 
stripping event. Upper panel: gas distribution in arbitrary 
units. The relative contours are the same as in Fig. |3 
Middle panel: velocity field. The contours are the same as 
in Fig. 13 Lower panel: velocity dispersion. The contours 
are the same as in Fig. O 
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Fig. 15. Three dimensional view of the data cubes. The 
RA-Dec plane is perpendicular to the image plane, and 
the velocity axis is the vertical axis, a) H i observations, b) 
Unperturbed gas disk, c) Ongoing ram pressure stripping, 
d) Past ram pressure stripping event. 



